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We study the heat kernel for the Laplace type partial differential op- 
erator acting on smooth sections of a complex spin-tensor bundle over a 
generic n-dimensional Riemannian manifold. Assuming that the curvature 
of the U(l) connection (that we call the electromagnetic field) is constant 
we compute the first two coefficients of the non-perturbative asymptotic 
expansion of the heat kernel which are of zero and the first order in Rie- 
mannian curvature and of arbitrary order in the electromagnetic field. We 
apply these results to the study of the effective action in non-perturbative 
electrodynamics in four dimensions and derive a generalization of the 
Schwinger's result for the creation of scalar and spinor particles in elec- 
tromagnetic field induced by the gravitational field. We discover a new 
infrared divergence in the imaginary part of the effective action due to the 
gravitational corrections, which seems to be a new physical effect. 
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1 Introduction 

The effective action is one of the most powerful tools in quantum field theory and 
quantum gravity (see [fT8l[T5l l9l [T0l[T3l ). The effective action is a functional of the 
background fields that encodes, in principle, all the information of quantum field 
theory. It determines the full one-point propagator and the full vertex functions 
and, hence, the whole S -matrix. Moreover, the variation of the effective action 
gives the effective equations for the background fields, which makes it possible 
to study the back-reaction of quantum processes on the classical background. In 
particular, the low energy effective action (or the effective potential) is the most 
appropriate tool for investigating the structure of the physical vacuum in quantum 
field theory. 

The effective action is expressed in terms of the propagators and the vertex 
functions. One of the most powerful methods to study the propagators is the heat 
kernel method, which was originally proposed by Fock [|T6ll and later generalized 
by Schwinger |[T%1l who also applied it to the calculation of the one-loop effective 
action in quantum electrodynamics. Finally, De Witt reformulated it in the geo- 
metrical language and applied it to the case of gravitational field (see his latest 
book (151). 

In particular Schwinger solved exactly the case of a constant electromagnetic 
field and derived an heat kernel integral representation for the effective action. He 
showed that the heat kernel becomes a meromorphic function and a careful evalu- 
ation of the integral leads to an imaginary part of the effective action. Schwinger 
computed the imaginary part of the effective action and showed that it describes 
the effect of creation of electron-positron pairs by the electric field. This effect 
is now called the Schwinger mechanism. This is an essentially non-perturbative 
effect (non-analytic in electric field) that vanishes exponentially for weak electric 
fields. 

Therefore its evaluation requires non-perturbative techniques for the calcula- 
tion of the heat kernel in the situation when curvatures (but not their derivatives) 
are large (low energy approximation). A powerful approach to the calculation of 
the low-energy heat kernel expansion was developed in non-Abelian gauge theo- 
ries and quantum gravity in|[21IS[5j[S|7l[M[IIl[l2]|. While the papers E H 
dealt with the constant electromagnetic field in flat space, the papers flU [7J O 
dealt with symmetric spaces (pure gravitational field in absence of an electromag- 
netic field). The difficulty of combining the gauge fields and gravity was finally 
overcome in the papers IfTTl [T211 . where homogeneous bundles with parallel cur- 
vature on symmetric spaces was studied. 
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In lfT4l we computed the heat kernel for the covariant Laplacian with a strong 
covariantly constant electromagnetic field in an arbitrary gravitational field. We 
evaluated the first three coefficients of the heat kernel asymptotic expansion in 
powers of Riemann curvature R but in all orders of the electromagnetic field F. 
This is equivalent to a partial summation in the heat kernel asymptotic expansion 
as t —> of all powers of F in terms which are linear and quadratic in Riemann cur- 
vature R. In the present paper we use those results to compute explicitly the terms 
linear in the Riemann curvature in the non-perturbative heat kernel expansion for 
the scalar and the spinor fields and compute their contribution to the imaginary 
part of the effective action. In other words, we generalize the Schwinger mech- 
anism to the case of a strong electromagnetic field in a gravitational field and 
compute the gravitational corrections to the original Schwinger result. 



2 Setup of the Problem 

Let M be a n-dimensional compact Riemannian manifold (with positive-definite 
metric g^ v ) without boundary and S be a complex spin-tensor vector bundle over 
M realizing a representation of the group Spin(n) ® t/(l). Let <p be a section of 
the bundle S and V be the total connection on the bundle S (including the spin 
connection as well as the U(l) -connection). Then the commutator of covariant 
derivatives defines the curvatures 

\y^y\<P = (^ + iF^)(p, (2.1) 

where F MV is the curvature of the f/(l) -connection (which will be also called the 
electromagnetic field) and 7? pv is the curvature of the spin connection defined by 

^fiv = /jv^ab ■> (2.2) 

with E flfo being the generators of the spin group Spin(n) satisfying the commutation 
relations 

[X ab X d ] =45 [c [a xV (2-3) 
Note that for the scalar fields ^ v = and for the spinor fields 

Z ab = ^Jab , (2.4) 

where y ab = y [a y b] (more generally, we define y ai ...a m = 7[a, ' ■ ' 7a m ]) and y a are the 
Dirac matrices generating the Clifford algebra 

JaJb + JbJa = 2g ab l . (2.5) 
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2.1 Differential Operators 

In the present paper we consider a second-order Laplace type partial differential 
operator, 

L = -A + £R + Q, (2.6) 

where A = g^ v V^V v is the Laplacian, £ is a constant parameter, and Q is a smooth 
endomorphism of the bundle S. This operator is elliptic and self-adjoint and has 
a positive-definite leading symbol. Usually, for scalar fields we set 

2 scalar = . (2.7) 

Moreover, for canonical scalar fields the coupling 

for canonical scalar fields , 

in - 2) ( 2 - 8 ) 
for conformal scalar fields . 



^scalar _ ^ 



I 4(n - 1) 

Another important case is the square of the Dirac operator acting on spinor fields 

L = D 2 , (2.9) 

where 

D = iy ti V M . (2.10) 
It is easy to see that in this case we have 

^spinor = 1 (2n) 

and 

Q spinor = -^/>/ v . (2.12) 



2.2 Effective Action 

The object of primary interest in quantum field theory is the (Euclidean) one-loop 
effective action determined by the formal determinant 

r ( i) = crlogDet(L + m 2 ), (2.13) 

where <x is the fermion number of the field equal to (+1) for boson fields and (-1) 
for fermion fields, m is a mass parameter, which is assumed to be sufficiently large 
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so that the operator (L + m 2 ) is positive. Notice that the usual factor \ is missing 
because the field is complex, which is equivalent to the contribution of two real 
fields. Of course, this formal expression is divergent. To rigorously define the 
determinant of a differential operator one needs to introduce some regularization 
and then to renormalize it. One of the best ways to do it is via the heat kernel 
method. 

In a physical theory the effective action describes the in-out vacuum transition 
amplitude via 

(out|in) = exp[/r (1) ] . (2.14) 

The real part of the effective action describes the polarization of vacuum of quan- 
tum fields by the background fields and the imaginary part describes the creation 
of particles. Namely, the probability of particles production (in the whole space- 
time) is given by 

P= 1 - Kout|in}| 2 = 1 -exp[-2Imr (1) ] . (2.15) 

Of course, the unitarity requires that the imaginary part of the effective action 
should be positive 

lmr (1) >0. (2.16) 

Notice that usually, when the imaginary part of the effective action is small, we 
just have 

/ , «2Imr ( i ) . (2.17) 
The one-loop effective Lagrangian is defined by 

r (1) = J dxg ll2 £. (2.18) 

M 

Therefore, the rate of the particles production per unit volume per unit time is 
given by the imaginary part of the effective Lagrangian 

R = — «2Im£. (2.19) 
VT 

2.3 Spectral Functions 

The heat kernel for the operator L is defined as the solution of the heat equation 



(d, + L) U(t\x, x') = 0, 



(2.20) 



/. G. Avramidi and G. Fucci : Non-perturbative QED in Curved Spacetime 5 



with the initial condition 

U(0\x,x') = 8{x,x') . (2.21) 

where 6(x, x') is the covariant scalar delta function on the bundle S. The heat 
kernel diagonal is defined by 

U diag (t) = U(t\x,x). (2.22) 

One of the best ways to describe the spectral properties of the operator L is via 
the heat trace 

Tr exp(-tL) = J dx g 1,2 Q(t) , (2.23) 

M 

where dx is the Lebesgue measure on the manifold M, g = detg^v and 

0(0 = tr f/ diag (0 . (2.24) 

Here tr denotes the fiber trace over the bundle S. 

The determinant of the operator can be defined within the so-called zeta- 
function regularization as follows. First, one defines the zeta function by 

£(s) = fi 2s Tr (L + m 2 \ S = J dx g U2 Z(s) , (2.25) 

M 

where 



Z(s) = ^--jdt rV tm2 0(O , (2.26) 
o 

and p is a renormalization parameter introduced to preserve dimensions. The zeta 
function £(s) is a meromorphic function of s analytic at s = 0. This enables one 
to define the (zeta-regularized) functional determinant of the operator (L + m 2 ) by 

Det (L + m 2 ) = exp[-^'(0)] , (2.27) 

where £'(s) = j-^(s). Therefore, the one-loop effective action is simply 

r (1) = -<(0), (2.28) 

and the one-loop effective Lagrangian is given by 

L = -o-Z'(O) . (2.29) 
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The effective Lagrangian can be also defined simply in the cut-off regularization 
by 




(2.30) 



where s is a regularization parameter, which should be set to zero after subtracting 
the divergent terms. Another regularization is the dimensional regularization, in 
which one simply defines the effective action by the formal integral 




(2.31) 



o 



where the heat trace is formally computed in complex dimension (n - 2s) with 
sufficiently large real part of s so that the integral is finite. The renormalized 
effective action is obtained then by subtracting the simple pole in s. 

For elliptic operators (in the Euclidean setup) the heat trace is a smooth func- 
tion of t; in many cases it is even an analytic function of t in the neighborhood 
of the positive real axis. However, in the physical case for hyperbolic operators 
(in the Lorentzian setup) the heat trace can have singularities even on the positive 
real axis of t. As we will show later in the approximation under consideration (for 
constant electromagnetic field) it becomes a meromorphic function of t with an 
essential singularity at t = and some poles f ft , k = 1, 2, . . . „ on the positive real 
axis. It turns out that the imaginary part of the effective action does not depend 
on the regularization method and is uniquely defined by the contribution of these 
poles. These poles should be avoided from above, which gives 

CO 

Im£ = -crnj^ Res [t~ l e~ tml ®{t)\ t k ) . (2.32) 

k=\ 

This method was first elaborated and used by Schwinger [18J in quantum electro- 
dynamics to calculate the electron-positron pair production by a constant electric 
field. One of the goal of our work is to generalize the Schwinger results for the 
case of constant electromagnetic field in a gravitational field. We will compute 
the extra contribution to the particle production by a constant electromagnetic 
field induced by the gravitational field. 
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2.4 Heat Kernel Asymptotic Expansion 

It is well known iTTVTl that the heat kernel diagonal has the asymptotic expansion 
as t -> (see also EEC SOS) 

CO 

U**(t) ~ {Anty nl2 Yj t k <*k , (2.33) 

k=0 

where a* are the local heat kernel coefficients. Then the trace of the diagonal heat 
kernel has the corresponding asymptotic expansion 

oo 

G(t)~(4ntr n,2 J]t k A k , (2.34) 

k=0 

where 

A k =tra k . (2.35) 

The diagonal heat kernel coefficients a, are polynomials in the jets of the met- 
ric, the U(l) -connection and the potential term Q, that is, in the curvature tensors, 
Q, and their derivatives. The lower order diagonal heat kernel coefficients are well 

known ei mm 

a Q = 1 , (2.36) 

a, = -Q + 11-Ar- (2-37) 

To avoid confusion we should stress that the normalization of the coefficients a k 
differs from the papers (21 El 13. 

In our previous paper lfT4l we studied the case of a parallel U(X) curvature 
(covariantly constant electromagnetic field), i.e. 

Vtf = . (2.38) 

In the present paper we will also assume that the potential term Q is covariantly 
constant too 

V iU 2 = 0. (2.39) 

By summing up all powers of F in the asymptotic expansion of the heat kernel 
diagonal we obtained a new (non-perturbative) asymptotic expansion 

oo 

U**(t) - {AntV 12 exp (-tQ) J(t) J] ^h(t) , (2.40) 

k=Q 



/. G. Avramidi and G. Fucci : Non-perturbative QED in Curved Spacetime 8 



where 

/ tiF \ m 
7(0 = det — -— (2.41) 
\ sinh(nF) / 

and bk(t) are the modified heat kernel coefficients which are analytic functions of 
t at t = which depend on F only in the dimensionless combination tF . Here and 
everywhere below all functions of the 2-form F are analytic at and should be 
understood in terms of a power series in the matrix F = (F^ v ). Notice the position 
of indices here, it is important! There is a difference here between Euclidean case 
and the Lorentzian one since the raising of indices by a Minkowski metric does 
change the properties of the matrix F . Also, here det denotes the determinant with 
respect to the tangent space indices. 

The fiber trace of the heat kernel diagonal has then the asymptotic expansion 

oo 

0(0 ~ (4nty n/2 ®(t) J] • ( 2 - 42 ) 

k=0 

where 

3>(O = /(OtrexpH0 , (2.43) 
tr exp (-tQ) 

are new ( non-perturbative) heat kernel coefficients of the operator L. The integrals 
f dxg l/2 B k (t) are then the spectral invariants of the operator L. 

This expansion can be described more rigorously as follows. We rescale the 
£/(l)-curvature F by 

F i — > F(t) = r l F , (2.45) 

so that tF(t) = F is independent of t. Then the operator L(t) becomes dependent 
on t (in a singular way!). However, the heat trace still has a nice asymptotic 
expansion (12.421) as t — > 0, where the coefficients are expressed in terms of 
F = tF(t), and, therefore, are independent of t. Thus, what we are doing is the 
asymptotic expansion of the heat trace for a particular case of a singular ( as 
? — > time -dependent operator L{t). 



3 Calculation of the Coefficient B\{t) 

For the first two coefficients we obtained ffT4l 

b (t) = 1, (3.1) 
b,(t) = {^^(0 + ^(0}^ (3.2) 
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where 



W(t) = I (coth(riF) - 1 , (3.3) 



2\ tiF 



i 



-^^[v(t)^](t) - l -^\{T)^ a \r)^ , (3.4) 

and 

1 exp[(l - 2r)riF] - exp(-riF) 

= o • rTTT^r , (3.5) 

coth(riF) 1 

S(t) = ~ - F ■ , . . F , cosh[(l-2T)riF], (3.6) 
tiF tiF smh(tiF) 

Z(t) = 3ft"Fcoth(ft'F) + ^T" — cosh[(l - 2r)ft'F] . (3.7) 

sinn(/iF) 

The trace coefficients are then given by 

B (t) = 1, (3.8) 

Bi(t) = (V)^(0 + V^(f)}^, (3-9) 



where 

tr exp (-J0 
tr exp i-tQ) 



vau = —^-7^7^- (3.io) 



3.1 Spectral Decomposition 

To evaluate it we use the spectral decomposition of the matrix F = (F M V ), 

N 

F = Y J B k E k , (3.11) 

k=\ 

where B k are some real invariants and E k = (E k M v ) are some matrices satisfying 
the equations 

E k fi V = ~E k yfj , (3.12) 
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£^=0, (3.13) 

and fork^m 

E k E m = 0. (3.14) 

Here, of course, N < [n/2]. The invariants B k (that we call "magnetic fields") 
should not be confused with the heat trace coefficients 5 and B\. 
Next, we define the matrices 11* = (11/ v ) by 

U k = -E 2 k . (3.15) 

They satisfy the equations 

n fejUV = n feVAt , (3.16) 

E k U k = U k E k = E k , (3.17) 

and forkj^m 

E k U m = U m E k = , U k U m = . (3.18) 

To compute functions of the matrix F we need to know its eigenvalues. We 
distinguish two different cases. 



Euclidean Case. In this case the metric has Euclidean signature (+ + h) and 

the non-zero eigenvalues of the matrix F are +iB u +iB N , (which are all imag- 
inary). Of course, it may also have a number of zero eigenvalues. In this case the 
matrices n* are nothing but the projections on 2-dimensional eigenspaces satisfy- 
ing 

(3.19) 
(3.20) 



n z k = u k 



In this case we also have 



B t = -E^ V F, 



k k lv ' 



Then we have 



N 

(iF) 2m = Y J B 2 k m U k , (m>l) 

k=\ 

N 

2m+l 



Y,B 2 rHE k , (m>0), 



and, therefore, for any analytic function of tiF at t = we have 

N 



k=\ 



f(tiF) = /(0)I + £ ~\f(tB k ) + f{-tB k ) - 2/(0) 



1 

IT + - 



f(tB k )-f(-tB k ) 



(3.21) 
(3.22) 

(3.23) 
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Pseudo-Euclidean Case. This is the physically relevant case of pseudo-Eucli- 
dean (Lorentzian) metric with the signature (- + h). Then the non-zero eigen- 
values of the matrix F are ±B\ (which are real) and ±iB 2 , . . . , ±iB N , (which are 
imaginary). We will call the invariant B\, determining the real eigenvalue, the 
"electric field" and denote it by B\ = E, and the invariants B k , k = 2, . . . , N, de- 
termining the imaginary eigenvalues, the "magnetic fields". So, in general, there 
is one electric field and (N - I) magnetic fields. Again, there may be some zero 
eigenvalues as well. 

In this case the matrices n 2 ,. . . JIn are the orthogonal eigen-projections as 
before, but the matrix IT is equal to the negative of the corresponding projection, 
in particular, 



n? = -n l5 u.e^-e,, n^ = -2. 



Now, we have 



N 



k=2 
N 



(3.24) 



(3.25) 



(iF) 2m+l = (iE) 2m+l E, + J] B 2 k m+l iE k , (m > 0) , (3.26) 



k=2 



Thus, to obtain the results for the pseudo-Euclidean case from the result for 
the Euclidean case we should just substitute formally 



B\ i-> iE, iE\ i-> Ei, III i-^ ~rii 
In this way, we obtain for an analytic function of itF, 



(3.27) 



f(tiF) = f (0)1- ^f(itE) + f(-itE) -2/(0) 
+ J]\-[f(tB k ) + f(-tB k ) -2f(0) 



1 

n, + - 



l 

n, + - 



fdtE)-f(-itE) 
f(tB k )-f(-tB k ) 



iE k 
(3.28) 



3.2 Scalar and Spinor Fields 

First of all, we note that for scalar fields 



o scalar (0 = j(t) . 



yscalar^ = q 



(3.29) 
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For the spinor fields we have 

*-« W = ,(-)Te X p(I„w) (3.30) 

r^J^M. (3.31) 
" 2 trexp(i^F pcr yP-) 

Here tr denotes the trace with respect to the spinor indices. 

We will compute these functions as follows. We define the matrices 

T k = l -iE^ v . (3.32) 
Then by using the properties of the matrices E k and the product of the matrices 

/ v r^ = - 4<Sjjr v] fl - 26%6 v fl , (3.33) 

(and some other properties of Dirac matrices in n dimensions) one can show that 
these matrices are mutually commuting involutions, that is, 

Tl = l, (3.34) 

and 

[T k ,T m ] = 0. (3.35) 
Also, the product of two different matrices is (for k ± m) 

T k T m = ~E^ v E^ 7Mvap . (3.36) 
More generally, the product of m > 1 different matrices is 

T kl --- T km = (-) E^ ■ ■ ■ E^y Mm . (3.37) 

It is well known that the matrices y m ...^ k are traceless for any k and the trace of 
the product of two matrices y Ml ... Mk and y Vl ... Vm is non-zero only for k = m. By using 
these properties we obtain the traces 

trT^O, (3.38) 

tryrfTk = -2 {n,1 \Ef , (3.39) 
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and for m > 1 : 

tr TV- -7^=0, (3.40) 
7^=0. (3.41) 

when all indices k\ , . . . k m are different. 

Now, by using the spectral decomposition of the matrix F we easily obtain 
first 

N _ 

11 smh(tB k ) 

and 

II \ N 

tr exp ^-tiF^j = tr II exp (tTkBk) ■ (3 - 43) 

/ 1 \ * 

tr^ exp \^tiF, v -f v j = tr fl ex P (tTkBk) ■ (3M) 

By using the properties of the matrices T k we get 

exp (tT k B k ) = cosh(^) + T k sinh(^) . (3 .45) 

Therefore 

/ i \ N 

tr exp {^tiF^J = 2 [n/2] [J cosh(^) , (3.46) 

and 

N N N 

try^Y\^{tT k B k ) = [J cosh(tBj) £ Umh(tB^try^T k 

k=\ j=l k=i 

N N 

= -2 [n/2] Y\ cosh(tBj) ^ tmh(tB k )iEf . (3.47) 

7=1 k=l 

Thus for the spinor fields 

N 

O spinor (0 = 2 [n/2] ]~[ ffl* cothftfl *) , (3.48) 



k=i 



and 

1 W 

^7° r (') = - o Z tanh (^)^^ • ( 3 - 49 ) 



k=i 
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By the way, this simply means that 

v F spinor (0 = ~ tanh(ft'F) . (3 .50) 

3.3 Calculation of the Tensor V^ av p(t) 

Next, we compute the tensor V^ypit). First, we rewrite in the form 

i 

o 

+1^(t)^(t) - lj^ [v ( T )jr ] ffl (T)| , (3.51) 

where 

*( T ) = -co t h(„-F) + COSh[(1 - 2 I ) " F1 , (3.52) 

sinh(^F) 

sinh[(l - 2r)tB k ] 

J/(r) = 1+ L . r - - (3.53) 

sinh(m Jt ) 

Next, we parametrize these matrices as follows 

AT 

S(t) = 2r(l-r)I + 2/u(T)n,, (3.54) 

yt=l 

/V 

Z(t) = 4I + 2/2jt(T)n*, (3.55) 

AT 

J/(r) = 2(l-r)I + ^A.(r)n,, (3.56) 
#(t) = ^/ 4 «, (3.57) 

N 

W(t) = J]f 5 , k (t)iE k , (3.58) 



4=1 
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where 



coth(tB k ) 1 
/u(r) = ' - - - . r - =rr cosh[(l-2r)^]-2r(l-r),(3.59) 

ffi^ sinh(^) 

/ u (t) = 3^coth(ffi0 + . ^- - cosh[(l - 2r)ffl fc ] - 4 , (3.60) 
sinh[(l - 2r)tB k ] 

hk( T ) = ■ 7— p. : (1-2t), (3.61) 

sinh(?6^) 

/ , ^ ^ ^ COsh[(1 ~ 2t) ^ ] <i ^ 

f 4 , k (r) = -coth(tB k ) + . , (3.62) 

sinh(?5j.) 

/ 5 ,,(0 = lLth(^)-^-j . (3.63) 



This parametrization is convenient because all functions f m , k (T) are analytic func- 
tions of t at t = and f m jjj) 
Then we obtain 



= 0. 

r=0 



iV N 

+ Z Z [p^(0n^ [v n| m | a % - o- km wrE m vp ] , (3.64) 



k=l m=l 

where 



<Pk(t) = -J2j dr i 2 f^(r) + r(l - T)f u (r) + 4(1 - r)/ 3) ,(r)] , (3.65) 
o 

l 

J rfT[/ 1 , fe (T)/ 2 , m (T)+/ 2 ,,(T)/ 1 , m (T) + 4/ 3) ,(T)/3 >m (T)] (3.66) 





1 



f*m(f) = J£ J dTf 4ik (T)f 4tm (T) . (3.67) 
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3.4 Calculation of the Coefficient Functions 

The remaining coefficient functions ifk(t), Pkmif) and cr km (t) are analytic functions 
of t at t = 0. To compute them we use the following integrals 

i 

/sinhx 
Jrcosh[(l -2t)x] = , (3.68) 

o 

i 

j Jrsinh[(l -2t)jc] = 0. (3.69) 
o 

By differentiating these integrals with respect to x we obtain all other integrals we 
need 

i 

/I sinh x 
drr cosh[(l -2t)x] = , (3.70) 
2 x 

o 
i 

/dr r 2 cosh[(l - 2t)x] = - - + — r sinhx - -—? coshx, (3.71) 
2\x x i j 2x l 



o 

l 



/. , „ , , lcoshjc 1 sum* 

rfrrsinh[(l-2T)x] = --— — + -— — . (3.72) 

o 



We also have the integrals 



i 

Jrcosh[(l - 2t)x] cosh[(l - 2r)y] = -< — + — } , 

2 [ x + y x-y j 

o 

(3.73) 

i 

J~dr cosh[(l - 2t)*] sinh[(l - 2r)y] = 0, (3.74) 



o 

l 



f a • urn o i • urn o ^ i 1 f sinh(* + y) sinh(x-y)\ 

dr sinh[(l - 2t)x] sinh[(l - 2r)y] = -< > 

J ' 2{ x+y x-y 



(3.15) 
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By using these integrals we obtain 

<Pk(t) 



13 11 / 1 

- + - - — coth(^) \tB k + 9 



CTkmif) = 



6 8 (tB k ) 2 24 



1 



— coth(^)coth(?5 m ) 
16 16 (o, 



(tB k ), 

1 coth(tB k ) 1 coth(tB m ) 



(3.76) 



16 tBv 



1 coth(?5 m ) + coth(^) 1 coth(?fi m ) - coth(^) 



32 tffl* + fi m ) 



32 



(3.77) 



Pkm(t) - 



_l/ 4 + 9 ^- + 9^- . 

48 \ (tf t ) 2 (tB m ) 2 tB k 



8-^- coth(?5 (t ) - 8 — coth(?5 m ) 



5. 



-(?£*) coth(^) - (tB m ) coth(tB m ) coth(tB k ) 



B„ 



B k B m 



-3-^ coth(tB m ) + 3 | + |coth(?5 m ) coth(^) 



B m | B k 
B k B m 

B k B m 
— + — +4 

B m B k 



B m B k 



coth(tB m ) + coth(tB k ) 



t(B k + B m ) 
coth(tB m ) - coth(tB k ) 



t(B k - B m ) 



(3.78) 



3.5 Trace of the Heat Kernel Diagonal 

The trace of the heat kernel diagonal in the general case within the considered 
approximation is given by 

0(0 ~ (4ntT n/2 (b(t) {1 + tB x {t) + •••}, (3.79) 

where the function O(0 was computed above and the coefficient Bi is given by 



Bi(t) = I ^ - Ar + £ {^ a (t)f 5 , k (t)iE v k % av p + ^(onf *py} 

' ' k=l 

N N 

+ Z Z \Pkmit)K n-Aavfl ~ aUmT^t^) ■ (3.80) 



k=\ m=l 

Let us specify it for the two cases of interest 
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3.5.1 Scalar Fields 

For scalar fields we have 



N 



o scalar (0 = j ] 



(3.81) 



k J sinh(^) 

/ i \ N 

B scalar (0 = [-- A R + £ <p k {t)lf*R„ v (3.82) 
' ' k=l 

N N 



k=l m=l 



3.5.2 Spinor Fields 

For the spinor fields we obtain 



N 



O spinor (0 = 2 ln/2] Y\tB k coth(tB k ), 



k=\ 



12 k=i 



(3.83) 



(3.84) 



N N 



k=l m=l 



where 



1 tanh(?5 m ) 1 tanh(^) 

hmit) = <r k Jt) + - — — + o— ^ 

o trfk o tD m 

-\ tanh(?5 m ) coth(^) - \ tanh(^) coth(*B m ) . (3.85) 
8 8 

3.6 Equal Magnetic Fields 

We will specify the obtained result for the case when all magnetic invariants are 
equal to each other, that is, 



B\ — • • • — Bfj — B 



(3.86) 
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3.6.1 Scalar Fields 

For scalar fields it takes the form 



\ sinh(?5) / 

+p{t)H^ H a f R^yp - o-(t)X^ v X^R Mav p , (3.88) 



where 



k=l k=l 

13 1 1 , 3coth(?5) 

<*> = S^S-^llT' (3 - 90) 

16 32 r5 32 S inh 2 (?5) 

5 3 1 1 „ , , _ 7 coth(Jfi) 1 1 ,„ _ N 

p(t) = + —tB coth(tB) + — (3.92) 

^' 24 S(tB) 2 24 V 16 tB 16 sinh^fl) 

3.6.2 Spinor Fields 

For the spinor fields we obtain 

O spinor (0 = 2 M2] [tBcoth(tB)f , (3.93) 

5 s r or (o = 

+p(t)H i ; v HfR fiav/3 - A(t)X7*T R ^p , (3.94) 



where 



3 3 1 1 tanh(?5) 3 coth(?5) „ _ 

*° = + (3 - 95) 
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3.7 Electric and Magnetic Fields 

Now we specify the above results for the pseudo-Euclidean case when there is 
one electric field and (N - 1) equal magnetic fields. By using the recipe (13.271 ) we 
obtain the following results. 



3.7.1 Scalar Fields 

For scalar fields we have 

tE 



o scalar (o 



tB 



sm(tE) \sinh<>5) 



N-l 



(3.96) 



5f alar (0 



+a(t)E tl "E v { i R liav i } - Ip^H^IffR^avp + 2<x i(t)X™E v fR Mav p 



+p(t)H^' 'H^ 'R^p - cr{t)X'^ t X v 2 R^p 



(3.97) 



where 



N 



N 



k=2 



k=2 



(3.98) 



Pit) 

o-i (r) 



13 1 1,^3 cot(f£) 

tE cot(tE) + — — , 

6 S(tE) 2 24 v 8 tE 

5 3 1 1 , ^ 7 cot(tE) 1 

-— + -t=tt + —tEcot(tE) - — + 



(3.99) 



1 



24 8 (tE) 2 24 



1 3 cot(rE) 3 



16 16 sin 2 (tE) ' 

(3.100) 



1 



16 32 tE 



32 sin 2 (?£) ' 

1 cot(tE) 1 coth(ffl) 



— cot(f£) cofh(?5) 

16 v v 16 tB 16 r£ 



(3.101) 



+- 



1 B cot(tE) + E coth(ffl) 



16 



?(5 2 + £ 2 ) 



(3.102) 
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Pi(0 = 



--U - 9— !-r + 9-^ + 8— cot(tE) - 8— coth(Jfi) 
48 \ (tE) 2 (tB) 2 tE tB 



-(tE) cot(tE) - (tB) coth(*£) 
E . _ B 



-3 — cot(^) + 3— coth(?5) + 3 (- - - J coth(?5) cot(tE) 



5B 2 -E 2 , , „ N 
+^r-^ — — coth(ffl) 



?5(5 2 + £ 2 ) 



5£ 2 - fl 2 
?£(fl 2 + E 2 ) 



cot(tE) \ . 



(3.103) 



3.7.2 Spinor Fields 

For the spinor fields we obtain 

<D spinor (0 = 2 [n/2] tEcot(tE)[tBcoth(tB)f- [ , 



(3.104) 



B s r°\t) = ~R - y(t)K R ^ + <PW?R MV + p(onf nf 

+A(t)El a E y fR tia yi 3 - Ip^fjH^TffR^avp + 2Ai(t)X^ a E v f R^p 
+p{t)H l ^H°fR liav i 3 - Ai^X^X^R/jayp , 



where 



Xt) 
Ut) 



1 tan(jE) 3 cot(*E) 



(3.105) 
(3.106) 



16 32 srn 2 (;£) 4 ?E 32 ' 

— cot(*E) coth(ffi) - — — — — — — 

16 16 tB 16 tE 

1 Bcot(tE) + Ecoth(tB) 1 tanh(Jfi) 1 tan(jE) 

+ 16 ?(fl 2 + £ 2 ) + 8 tE 8 tB 

1 1 

-- tanh(ffi) cot(*E) + - tan(*E) coth(?fi) . (3. 107) 

8 8 



4 Imaginary Part of the Effective Lagrangian 

Now, we can compute the imaginary part of the effective Lagrangian in the same 
approximation taking into account linear terms in the curvature. The effective ac- 
tion is given by the integral over t of the trace of the heat kernel diagonal. Of 
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course, it should be properly regularized as discussed above. The most important 
point we want to make is that in the presence of the electric field the heat kernel is 
no longer a nice analytic function of t but it becomes a meromorphic function of t 
in the complex plane of t with poles on the real axis determined by the trigonomet- 
ric functions in the coefficient functions computed above. As was pointed out first 
by Schwinger these poles should be carefully avoided by deforming the contour of 
integration which leads to an imaginary part of the effective action determined by 
the contribution of the residues of the poles. This imaginary part is always finite 
and does not depend on the regularization. We compute below the imaginary part 
of the effective Lagrangian for the scalar and the spinor fields. 

The trace of the heat kernel 0(f) was computed above and is given by (13.791) . 
Now, by using (12.321 ) the calculation of the imaginary part of the effective La- 
grangian is reduced to the calculation of the residues of the functions t~ n l 2 ~ l e~ m O(f) 
and t~ n/2 e~ tm <J>(f).Bi(f) at the poles on the real line. By using the result (13.961 ) and 
(13.1041) for the function O it is not difficult to see that the function t~ nl2 ~ l e~ tn? Q>{t) 
is a meromorphic function with isolated simple poles at f# = kn/E with k = 

1,2 The function r n/2 e~ tm O(Y)5i(0 is a ls° a meromorphic function with the 

same poles but the poles could be double or even triple. The imaginary part is, 
then, simply evaluated by summing the residues of the integrand at the poles. It 
has the following form 



Im£ = n(4nr n/2 E n/2 G (x,y) + n(4ny n,2 E n/2 - 1 

+ G 2 (x, y)IIf + G 3 (x, y)H/ v R MV + G 4 (x, y)nf If;"/? 

<1 f ^S{X, y)ni n 2 ^Rfiav/3 



Gi(x,y)R 

fiavp 

+ G 5 (x, y)E f ; a EfR IJavf3 + G 6 (x, yW v Iff R, 
+ G 7 (x, y)X 2 » a EfR, avl3 + G 8 (jc, y)/// v // 2 ^ 



where 



+ G 9 (x,y)X 2 ^ a X 2 v % av ^ 

B m 2 
E ' y = l 
and Gj(x, y) are some functions computed below 



(4.1) 



x=-, y = T , (4.2) 



4.1 Scalar Fields 



At this point it is useful to introduce some auxiliary functions so that the final 
result for the quantities G^ calm (x, y) can be written in a somewhat compact form, 
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namely 



fk(x,y) = 



knx 



N-l 



exp (-kny) 



sinh (knx) 

gk(x,y) = (N - \)(knx) coth(knx) + kny , 

h k (x, y) = ^N(N - l)(knx) 2 coth 2 (fon:) + ^ - /v) kny 

+ ^(N - l)[(n - 2N) + 2kny](knx) coth(for;t) 
1 



(4.3) 
(4.4) 



+ -(kn) 2 [l-(N-l)x 2 +y 2 ] 



(4.5) 



k(x, y) = -knx + 



(i ~ N ) + kny 



coth(for;t) + N(knx) coth 2 (fou;) , (4.6) 



1 n-2N 



1 



S(kn) 2 1 2 / 



+— 1 - 
24 



(knf 



8k(x,y) , 



&2,k(x,y) = --- 



1 n-2N 
6 



48 



l 



21 \ 1 ^fcC*,?) 



2{kn) 



16 (for) 2 



""^ " T-6-6^{l- N+1 ){h N + 2 ) 



[h k (x,y) + g k (x,y)] , 



Q 4 , k (x,y) = --- 



32(kn) 2 
1 n-2N 



8 



48 



+ 



8(for) 2 



**(*,y) 

kn 



+ 



8(fc 



1 coth(fco:) fx 4 — 3 (for*) 2 



3* 4 - 1 



8 knx 
1 



+ 1 



8(jb#lV2- W+1 
1 coth(£:;o;) / 



1 + g*(*.;y) 



Jt 2 (;t 2 + 1) 



l 



1 l k (x,y) 
x(x 2 + 1) 8 (for) 



8 knx \ x 2 + 1 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 
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By using these quantities we obtain the functions G^ calar (x, y) in the form of 
the following series 



(—] ) k+l 

G^(x,y) = ^t-LpMx.y), (4.12) 

6-*jZhr^ Mx - yh (4 ' 13) 

G^(x,y) = -J^±—^f k ( x ,y)n uk (x,y), (4.14) 

G^(x,y) = J^±—^f k ( x ,yMk7rx), (4.15) 

k=\ 

G^(x,y) = Yj7}-^M x ^,k(x,y), (4.16) 

k=\ 

~ (-U k+1 

Gf^(x,y) = J^±-^ Mx ,y)n u ( x ,y), (4.17) 

(jC ,y) = _g^J_/ fc ( Xf y)n 4Jt ( Xf y), ( 4 .18) 

G^(x,y) = Y±-^—f k ( x ,y)n u (x,y), (4.19) 

G^Ocy) = ^ ±—^Mx,y)p(knx) , (4.20) 

G; calar (i,y) = V ; ? / /2 (4.21) 
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4.2 Spinor Fields 

Exactly as we did in the previous section, we introduce, now, some auxiliary func- 
tions that will be useful in the presentation of the final result, namely 

fs,k(x,y) = [(knx) coth(forx)] exp(-fory) , (4.22) 
8sjc(x, y) = (N- l)(knx) coth(knx) - (N - \){knx) tanh(knx) + kny , (4.23) 

h s , k {x,y) = \(kny) 2 -{N-\f{knxf + {^-Nyny (4.24) 

+ ^(N - 1) (n - 2N + 2kny) (knx)[ coth(knx) - tanh(fc^)] 

+ ^N(N - \)(knxf coih 2 (knx) + ^(N - l)(N - 2)(knx) 2 tanh 2 (for;t) , 

h,k{x, y) = -Nknx + ^-N + fory j coth(fc^) + N(knx) coth 2 (fc7rjt) , (4.25) 

Ps,k(x, y) = (N - 2)knx + ^-N + kjiyj tanh(ta:) - (N - 2)(knx) tanh 2 (to) , 

A , N 1 n-2N 3 In Ar \ 1 / 9 \ , N 



(4.26) 



, x 1 n-2N 1 In A7 \ln Ar \ 



(4.27) 



~ 24 l 1 -2wH JC,y) + T6lfa^-' (4 - 28) 

^ = -i|-6w(i-^ +1 )(i-^ +2 ) 

-p[gs,^,y) + /i w (x,y)] , (4.29) 



32(fcr) 
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+ 



1 n - IN 3 
8 48~~ + 

1 In 



S(kn) 2 



x 1 8 \ x 



kn 



3x 4 -l 



8(fc 

1 coth(fc^) 
8 fco; 



- 3 (fco;) 2 



+ 1 



x 2 (;c 2 + 1) 
1 



(4.30) 



— — - - N + 1 +g s , k (x,y)\ — — 

8(Aot) 2 \2 I x{x l + 1) 8 



kn 



1 cottsiknx) ( x 3 



8 knx \ x 2 



x 3 \ 1 

— I + -^[tmh(k7rx) + p s , k (x,y)] .(4.31) 



By using the above functions we can write the explicit expression for the quan- 
tities G* pmo W) 



oo 1 

Grw) = ^]g^ /w( ^), 



Gr°w) = - 



CO ^ 

Gf nor (x,y) = -2^J]—— I f s , k (x,y)A u (x,y), 

k=i ^ Kn) 

CO - 

or-two = 2 [51 E(^/»(-.>')^--). 

CO ^ 

co 1 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



(4.38) 
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DO ^ 

G s ^°\x,y) = 2^^^-^f s , k (x,y)A 5 , k (x,y), (4.39) 

oo 1 

en^) = 2[n Tj7j^^ fs ' k(x ' y)p(knx) ' (4 - 40) 

G; pi ™(^) = ~ 2[§] Z {kn) nl2^ k{x,yWknx) . (4.41) 

Notice that because of the infrared cutoff factor e~ kny the functions G ; (x, y) are 
exponentially small for massive fields in weak electric fields when the parameter 
is large, y » 1 (that is, m 2 » E), independently on x. In this case, all these 
functions are approximated by just the first term of the series corresponding to 
k=l. 

5 Strong Electric Field in Four Dimensions 

The formulas obtained in the previous section are very general and are valid in 
any dimensions. In this section we will present some particular cases of major 
interest. 

5.1 Four Dimensions 

In this section we will consider the physical case when n = 4. Obviously in four 
dimensions we only have two invariants, and, therefore, N = 2. The imaginary 
part of the effective Lagrangian reads now 

ImX = n(4ny 2 E 2 G (x, y) + n{Any 2 E G x (x, y)R 

+ G 2 (x, y)nf + G 3 (x, y)U/ v R, v + G 4 (x, y)tf? Iff R Mav/3 

+ G 5 (x, y)E t ; a EfR, av/s + G 6 (x, yW'Ilf R^ 

+ G 7 (x, y)ErE\%^ + G 8 (x, y)U/ v U 2 a % av/3 

+ G 9 (x,y)E 2 « a E 2 v % a J. (5.42) 
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For scalar fields in four dimensions the functions G* calar (x, y) take the form 

„-kny 



G^(x,y) = -V 

7T Z_i 



G\ cala \x,y) 



-^scalar 



(x,y) 



Gf aiar (x,y) 



G s 4 caiai (x,y) 



j-i k sinh(knx) 

1 \ y e^' 
6 *]*2jsinh(Jbrjc) ' 

CO 

-*y- 

A-/ si 



(5.43) 



(5.44) 



y -bry 



3 1 / 9 \ I 
i — — - + — far - — py + x coth(fco:)l > , 
^sinh(^)[8 4(for) 2 24 \ kn) u v J J 

(5.45) 



-toy 



1 3 1 
— I 



1 



J sinh(&7n:) (6 S(knx) 2 24 



- — fco:coth(fco;) - - 



3 coth(fco;) 



8 knx 



(5.46) 



-kny 



13 



+ 



13 



2 2 

x y 
+ — + 



^sinh(^)( 96 16(kn) 2 32 32 \16(for) 24 



&7T 



16 ■ 16far 24 *coth(ferx) + -coth(ferx)^ 



(5.47) 



Gr lar fey) = 



1 



^ sinh(to) (64 32(kn) 2 
3 / 1 \ . . 3 



3^ 
64 



64 



32kn 



- ^ I y + t- I coth(fc7T.x:) - — x 2 coth 2 (fco:) 
32 T for/ 32 



(5.48) 



+ 



3x 2 



+ 



k J sinh(forjc) { 4 2(knx) 2 A(kn) 2 8 
kn 1/1 \ 1 

s-* 2+1) -8s^" 3j: r3 <x4 " 1)c<>,h2<ta) 

l 



+ 7T(* 2 " 1) 



4&7T 



x/ 8 U 



COth(fco:) > , 



(5.49) 
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x 2 + i JE s inh(to) [S(kn) 2 ' 8 



+ — (x 2 + 1) + 



Skn 



y( l +x 2)_ 1(1-^) 



8 

^ 2 (x 2 + l)coth 2 (to)l , 



coth(toc) 



(5.50) 



Gr u \x,y) 



-fey 



7 3 1 



^ sinh(toH 48 8 (to) 2 16 
coth(toc) > , 



1 2 
- — coth (toe) 



f kn 7 
+ 1 —x + 



24 16to 



(5.51) 



3 3 
coth(toc) + — coth 2 (toc) 

f-Jsinh(to)l 32 32to v 32 v 

fc=i v 



(5.52) 



For spinor fields in four dimensions the functions G; Pmor (x, y) take the form 

(5.53) 



Gr ot (x,y) = £ £ I coth(to)e^ , 



G s r°Uj) = -^coth(to)<T^, 



(5.54) 



Gf nor (x,j) = -4xVcoth(to)^ ---f-j 

+ \kn - -^-j \y + xcoth(toc) - xtanh(toc)] j , (5.55) 



Gf mo \x,y) 



113 1 1 

Ax V cotWto)*? - ^ - + - — tocoth(to) 

\6 8 (to) 2 24 



3 coth(toc) 
8 knx 



(5.56) 
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GT no \x,y) 



00 ( 1 13 

Ax y coth(knx)e- kny l - - + — — 
^ \ 6 16(tor) 



16 + 32 



J / 21 \ x 2 2 
24 r" 2^ + T6 COth ^ 



x / _ 21 3y 
24 r 2£tt 



- y J [ coth(^) - tanh(^)]| , (5.57) 



GT n °\x,y) 



= 4x ^ coth(fco;)e 



3 3x 2 
+ 



k=\ 



16 32(for) 2 32 



3j 3x 2 , 2 ,, » 
- coth 1 (knx) 



3f 

64 32(fcr) 32 
3~f (&7r + y ) [ coth ( knx ) ~ tanh (^"0] } . 



(5.58) 



+ 



— ^ co\h{knx)e 

k=i 



8 4(for) 2 \ 3* 2 



8 V 24 \ kn 24 1 forx 2 



+ 



+ 



_L + i.^ |T ._ 



1 for* / 9 

+ + — \kn 

Sknx 24 24 1 kn 



COth(k7TX) 

tanh(knx) 



- ^(x 4 - l)coth 2 (knx)\ , 



(5.59) 



A 00 ( 1 2 

crW) = _g coth(ta)e -.,|_ + L (x2 + 1) 



+ — - — (x 2 + 1) coth 2 (te) 
Skn 4 



x 

+ 8 



+ 8 



i-(l-.r)-y(l+.i 2 ) 

K7T 



— (1 +2x 2 ) + 2y{\ +x 2 ) 
kn 



coth(knx) 
tanh(knx) } , 



(5.60) 
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00 ( 7 3 1 1 

Gf D0 \x,y) = Ax y coth(forx)<rH - — - -- — rr + —knxcoth(knx) 
8 [48 S(knx) 2 24 

7 coth(forx) 1 , 2/I N ) 
+ T6^^-16 COtl,<ta) }' <5 ' 61) 

Gf"(x,y) = -4* £ coth(te) e -"-"{ - ^ + ^ coth 2 (ta) + 

I-— 1 v 



£=1 

3 coth(fc^) 
32 A:7rx 



(5.62) 



5.2 Supercritical Electric Field 

As we already mentioned above, the functions G^x, y) are exponentially small for 
massive fields in weak electric fields for large y = m 2 /E, as y — > 00. Now we are 
considering the opposite case of light (or massless) fields in strong (supercritical) 
electric fields, when y — > with a fixed x. This corresponds to the regime 

m 2 «B,E. (5.63) 



5.2.1 Scalar Fields 

The infrared (massless) limit for scalar fields is regular — there are no infrared 
divergences. This is due to the presence of the hyperbolic sine sinh(fco:) in the 
denominator, which gives a cut-off for large k in the series, and therefore, assures 
its convergence. The result for the massless limit in the scalar case can be simply 
obtained by setting y = in the above formulas for the functions G,(x,y). 



5.2.2 Spinor Fields 

The spinor case is quite different. The presence of the hyperbolic cotangent 
coth(fco;) does not provide a cut-off for the convergence of the series as k — > 00. 
This leads, in the spinor case in four dimensions, to the presence of infrared di- 
vergences as y = m 2 /E — > 0. By carefully studying the behavior of the series as 
k — > 00 for a finite y and then letting y — > we compute the asymptotic expansion 
of the functions Gf mo \x, y) as y — > 0. 
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We obtain 

= y + 0(y). (5-64) 

G s r\x,y) = -J-f + f + OCy), (5.65) 
3^ j 8 

GT D °\x,y) = -l* + ^ + 0(y), (5.66) 
3ny 4 

^spinor, , 1 J 2 ^ 2 x 3 xVx-24) + 18 

G 3 P (xj) = -- — - + - — + — log(7ry) + — + 0(y) , 

ony z 3ny 2n 12x 

(5.67) 

GT n °\x,y) = -^ + ^ + 0(y), (5.68) 
ony 8 

GfUy) = + ^ + (5-69) 

jr 2 — 1 r 2 7 X X 4 - 3 

67r(jr + l)y / 3^ y 2n(x z + 1) 

6n(9x 4 + 3x 2 - 4) - 36x(x 4 - 1) + tt 2 xV - 1) _ 
72^ + 1) + O0 °' 

= -^^lo g (^) + ^l^ + 0(,), (5.71) 
2n(x z + 1) \2n(x l + 1) 

^spinor, , 1 X 2 5 X 7 ^(jT* - 30) + 1 8 

G « ^ = ^-^-4^ l08(;ry) 72* + 0(y) > 

(5.72) 

Gr° r (x,);) = + f log(^y) - ^ + . (5.73) 

Thus, we clearly see the infrared divergences of order x 2 /y 2 = B 2 /m 4 , x/y = 
B/m 2 and log y = \og(m 2 /E). 



5.3 Pure Electric Field 



We analyze now the case of pure electric field without a magnetic field, that is, 
5 = 0, which corresponds to the limit x — > with fixed y. This corresponds to the 
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physical regime when 

B«m 2 ,E. (5.74) 

In this discussion we present the results in arbitrary dimension first and then we 
specialize them to the physical dimension n = 4. 

5.3.1 Scalar Fields 

We now evaluate the functions Gi(x, y) for x = and a finite y. In this limit we are 
presented with series of the following general form 

(--\\k+l p -kny 

xT u (y) = I ( " \ nl2 • (5-75) 
k=\ K 

This series can be expressed in terms of the polylogarithmic function defined by 

Li ; (z) = Z I ' (5 - 76) 

so that, we have 



xT lw (y) = -Un(-e^) . (5.77) 
It is not difficult to notice that the limit as x -> of the functions G s 3 calar , G* calar , 

G scalar 5 G scalar and (^scalar yanish identically, that is, 

G* calar (0, y) = G* calar (0, y) = G^(0,y) = G* calar (0, y) = G^ calar (0, y) = . (5.78) 

The explicit expression for the remaining non-vanishing G* calar for pure electric 
field in n dimensions is 

G s ^(0,y) = -n- n/2 Lh(-e- ny ) , (5.79) 

0.7) = -(i-fj^U,.^^), (5.80) 
Gl c ^(0,y) = -^i 7 ^{ 2 7r 3 yLi f _ 2 (- e ^) + (n + 4VLi | _ 1 (-e- 7r >') 

- 18^yLi f (- e -^)-9(n + 2)Li f+1 (-e^')}, (5.81) 
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Gf'CO.y) = 3 847 )»/2 + i { ~ ^yUn_ 2 (-e-^) - 4n 2 (2n + 9 - 3y 2 )Ii.- 1 (-g-^) 
+ 12(n + ^/ryU^-e^') + 3(n + 2)(n + 24)Li« +1 (-e^')} , (5.82) 

Gr lar (0,y) = 256 ^ /2+1 { 4 ^a ~ 3y 2 )Li f _!(- g -^) - 12n^Li f (- g -^) 

- 3n(n + 2)Li» + i(-e _;r> ')} . (5.83) 

In the physical case of n = 4 some of the polylogarithmic functions can be 
expressed in terms of elementary functions. In this case we have 

GT lar (0,y) = —Lhi-e-xy). (5.84) 
n 

G\ c ^(0,y) = -|I-^Iln(l+ e ->), (5.85) 

Gf-'CO.y) = ^(|^ + 8, 2 ln(l + e--) 

+ 18^vLi 2 (-e- 7rv ) + 54Li 3 (-e- 7rv )J , (5.86) 
1 r 1 6?i 3 vp~ ny 

G4 (0,y) = _{ T -Z_ + 4*2(17 -3v 2 )ln(l +e -) 

+ l92jzyLi 2 (-e- ny ) + 504Li 3 (-e" 7rv )} , (5.87) 

256* 3 I 

+ 72Li 3 (-e^')J . (5.88) 

We study now the behavior of these functions as y —> 0, which corresponds to 
the limit 

5 = 0, m 2 «E. (5.89) 
By taking the limit as y —> of the expression (|5.77l) and by noticing that 

Li n (-l) = -(l-2 1 -")^(n) , (5.90) 
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where £(x) denotes the Riemann zeta function, we obtain 

G^(0,0)= K nn!1 • (5-91) 

Next, by taking the limit as y — > and by using the formula (15.901) , it is not 
difficult to obtain 

G scalar (a()) = _ / * _ f \ _ jl-n^ {± _ j j ^ ( g q 2) 

^"(0,0) = _^i^{-(„ + 4)^(1 -2 2 "^g-l) 

+ 9(n + 2)(l-2-' !/2 )^ + l)}, (5.93) 



Gr' ar (0,0) = 38 J» /2+1 {^\2n + 9)(1 - 2 2 -" /2 ^(^ - l) 



- 3(n + 2)(n + 24)(l-2-" /2 )^ + l)}, (5.94) 

^(0,0) = ^{-^(l-^g-l) 

+ 3n(n + 2)(l-2- /2 K^ + l)}. (5.95) 

We consider, at this point, the physical case of four dimensions. By setting 
n = 4 in (15.911 ) we obtain 

G^ alar (0,0) = J_. (5.96) 
Now, we notice the following relation 

(l-2 2 -"/ 2 x(^-l) = ^(^-l), (5.97) 

where n(x) is the Dirichlet eta function. In the particular case of four dimensions 
we have that 

lim(l - 2 2 ^ n/2 )d^ -l) = 77O) = In 2 . (5.98) 

n->4 \ 2 / 
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By using the last remark we obtain the values of the functions G ( (n, y) in four 
dimensions 

G s 1 calar (0,0) = -|i-£jiln2, (5.99) 
G 2 calar (0,0) = JLln2- J^£(3), (5.100) 

0,0) = -^l n2 + ^(3)- (5-102) 
5.3.2 Spinor Fields 

For spinor fields the expressions for the non-vanishing G* pinor in the limit x — > 
are 

G ( f nor (0,v) = 2 [ " /2] ^" /2 Li, ^/2 (e- ^^> •) (5.103) 

2 [n/2] i 

G s r r (0^) = -— ^irLif-i(^) . (5-104) 

G 2 Pm ° r (0^) = - 5 g^{2^Li i _ 2( ^) + ( n + 4)^Li f _ l( ^) 

- 187ryLi»( e -^) - 9(n + 2)Li» +1 (e^')} , (5.105) 

G; pm ° r (0,j) = 384 ^/2 + i { ~ 16^0^2^-^) - 4^(2/1 + 12 - Sy 2 ^.-!^) 
+ 12(n + 12)7ryLi f (<?-*>') + 3(n + 2)(n + 24)Li f+1 (e- ?r:v )J , (5.106) 

Gf no \0,y) = ~2 [ " /2] 256 ^ w/2+1 {4^ 2 (4 + y^)U 8 -i( g ^) + 4«iyU.(^) 

+ n(n + 2)Li f+1 (e~ ?r> ')} . (5.107) 
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In the particular case of n = 4 the above results read 

Gf lor (0,)0 = ^U 2 (e-«>), (5.108) 

7T 

G spmor (() } = 1 H1 _ e -ny h ( g m) 

3n 

G! pmor (0,v) = L l fLIl -g^lnd-g^) 

- \%nyU 2 (e- ny ) - 5AU 3 {e' ny )^ , (5.110) 

Gr r (O,>0 = -^{^^-4^(20-3v 2 )ln(l-e--) 

- 192^yLi 2 (e- ?r> ') - 50403^')} , (5.111) 

Gf nor (0, y) = j|^p(4 + y 2 ) ln(l - e*>) - AnyU^') - 6Li 3 (^)} . 

(5.112) 



In the case of spinor fields, for n > 4, there is a well defined limit as y — » 0. In 
fact, by taking the massless limit, y — > 0, of the expression (15.1031) and noticing 
that 

Li„(l) = <T(n) , (5.113) 

we obtain 

G — (0 ,0) = ±^g) . (5.114) 

Analogously, in the limit as y — > the result for the remaining G* pmor can be 
written as follows 

2 [n/2] in \ 

Gf in ° r (0,0) = -±—n l - nl2 d^-l\, (5.115) 



12 12 

2 [«/2] 



<rw) = -4f^(<» +4 ^(M- 9( " +2)f (H}- 



(5.116) 
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+3(n + 2)(n + 24)<rg + l)}, (5.117) 

Gr r (0,0) = -2^ s J^{l6H» f g-l) + »0. + 2)fg + l)}. 

(5.118) 

We turn our attention, now, to the physical case of n = 4. From the expression 
in (|5.1 141) we obtain the following result 

Gr° r (0,0) = ^. (5.119) 

It is evident, from the expressions in (|5.109I )- (I5.1 121) , that the functions 
G^ pmoi (0,y) in four dimensions represent a special case since there is an infrared 
divergence as m — » (or y — > 0). This means that there is no well-defined value 
for the massless limit y — > 0. Instead, we find a logarithmic divergence, log(7ry). 
In order to analyze this case we set n = 4 from the beginning in the expressions for 
finite y, and then we examine the asymptotics as y — > 0. By using the equations 
(f57T09l )- (f57TT2T) we obtain 

G spinor (() } = lj og(ny) + 0(y)r (5.120) 

3n 

Gf no \0,y) = |-log(^y)-i- + ^(3) + OCy), (5.121) 
Gf nor (0,3;) = |-log(^)-i- + ^(3) + 0(y) > (5.122) 

Gf nor (0,v) = ^log(^) + ^(3) + 0(v). (5.123) 

Notice that, in four dimensions the functions G* pmor (x, y) are singular at the 
point x = y = 0. In particular, the limits x — > and y — > are not commutative, 
that is, the limits as x — » of the eqs. (|5.64l) - (15.73l) (obtained as y — » for a finite 
x) are different from the eqs. (I5.120I) - (I5.123I) (obtained as y — > for x = 0). 
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6 Concluding Remarks 

In this paper we have continued the study of the heat kernel and the effective action 
for complex (scalar and spinor) quantum fields in a strong constant electromag- 
netic field and a gravitational field initiated in IPT41 . We study here an essentially 
non-perturbative regime when the electromagnetic field is so strong that one has 
to take into account all its orders. In this situation the standard asymptotic ex- 
pansion of the heat kernel does not apply since the electromagnetic field can not 
be treated as a perturbation. In [14] we established the existence of a new non- 
perturbative asymptotic expansion of the heat kernel and computed explicitly the 
first three coefficients of this expansion. 

In the present paper we computed the first two coefficients (of zero and the first 
order in the Riemann curvature) explicitly in n-dimensions by using the spectral 
decomposition of the electromagnetic field tensor. We applied this result for the 
calculation of the effective action in the physical pseudo-Euclidean (Lorentzian) 
case and computed explicitly the imaginary part of the effective action both in the 
general case and in the cases of physical interest. We also computed the asymp- 
totics of the obtained results for supercritical electric fields. 

We have discovered a new infrared divergence in the imaginary part of the ef- 
fective action for massless spinor fields in four dimensions (or supercritical elec- 
tric field), which is induced purely by the gravitational corrections. This means 
physically that the creation of massless spinor particles (or massive particles in 
supercritical electric field) is magnified substantially by the presence of the gravi- 
tational field. Further analysis shows that a similar effect occurs for any massless 
fields (also scalar fields) in the second order in the Riemann curvature. This ef- 
fect could have important consequences for theories with spontaneous symmetry 
breakdown when the mass of charged particles is generated by a Higgs field. Such 
theories would exhibit a significant amount of created particles (in the massless 
limit an infinite amount) at the phase transition point when the symmetry is re- 
stored and the massive charged particles become massless. That is why this seems 
to be an interesting new physical effect that deserves further investigation. 

References 

[1] Avramidi I G 1990 The co variant technique for calculation of the heat kernel 
asymptotic expansion, Phys. Lett. B 238 92-97 



/. G. Avramidi and G. Fucci : Non-perturbative QED in Curved Spacetime 40 



[2] Avramidi I G 1991 The covariant technique for calculation of one-loop effec- 
tive action, Nuc. Phys. B 355 712-754. Erratum: Nucl. Phys. B 509 (1998) 
557-558. 

[3] Avramidi I G 1993 A new algebraic approach for calculating the heat kernel 
in gauge theories, Phys. Lett. B 305 27-34 

[4] Avramidi I G 1994 The heat kernel on symmetric spaces via integrating over 
the group of isometries, Phys. Lett. B 336 171-177 

[5] Avramidi I G 1995 Covariant algebraic calculation of the one-loop effec- 
tive potential in non-Abelian gauge theories and a new approach to stability 
problem, J. Math. Phys. 36 1557-1571 

[6] Avramidi I G 1995 Covariant algebraic method for calculation of the low- 
energy heat kernel, J. Math. Phys. 36 5055-5070. Erratum: J. Math. Phys. 
39(1998) 1720 

[7] Avramidi I G 1996 A new algebraic approach for calculating the heat kernel 
in quantum gravity, J. Math. Phys. 37 374-394 

[8] Avramidi I G 1999 Covariant techniques for computation of the heat kernel, 
Rev. Math. Phys., 11 947-980 

[9] Avramidi I G 2000 Heat Kernel and Quantum Gravity, (Berlin: Springer- 
Verlag) 

[10] Avramidi I G 2002 Heat kernel approach in quantum field theory, Nucl. Phys. 
Proc. Suppl, 104 3-32 

[11] Avramidi I G 2008 Heat kernel on homogeneous bundles, Int. J. Geom. 
Meth. Mod. Phys., 5 1-23 

[12] Avramidi I G 2008 Heat kernel on homogeneous bundles over sym- 
metric spaces, Comm. Math. Phys., DOI: 10.1007/s00220-008-0639-6, 
|arXiv:math/0701489[ 55 pp 

[13] Avramidi I G 2009 Mathemathical Tools for Calculation of the Effective 
Action in Quantum Gravity, in: Quantum Gravity, Ed. B. Booss-Bavnbek, 
G. Esposito and M. Lesch, (Berlin, Springer); arXiv:08 12.3363 [hep-th], 
71pp, 



/. G. Avramidi and G. Fucci : Non-perturbative QED in Curved Spacetime 41 



[14] Avramidi I G and Fucci G 2008 Nonperturbative heat kernel asymptotics on 
homogeneous Abelian bundles, arXiv:08 10.4889 [math-ph] 45 pp 

[15] De Witt B S 2003 The Global Approach to Quantum Field Theory. Oxford 
University Press, Oxford 

[16] Fock V A 1937 The proper time in classical and quantum mechanics, Izv. 
USSR Acad. Sci. Phys. 4-5, 551-568 

[17] Gilkey P B 1995 Invariance Theory, the Heat Equation and the Atiyah-Singer 
Index Theorem. (Boca Raton: CRC Press) 

[18] Schwinger J S 1951 On gauge invariance and vacuum polarization, Phys. 
Rev. 82, 664-679 

[19] Vassilevich D V 2003 Heat kernel expansion: user's manual, Phys. Rep. 388 
279-360 



